Accepted for Publication in the Astronomical Journal: March 25, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 

A LARGE SAMPLE OF BL LACS FROM SDSS AND FIRST 

Richard M. Plotkin,^ Scott F. Anderson/ Patrick B. Hall,^ Bruce Margon,'* Wolfgang Voces, ^ 
Donald P. Schneider,'' Gregory Stinson,'^''^ and Donald G. York' 
Accepted for Publication in the Astronomical Journal: March 25, 2008 

ABSTRACT 

We present a large sample of 501 radio-selected BL Lac candidates from the combination of the 
Sloan Digital Sky Survey (SDSS) Data Release 5 optical spectroscopy and from the Faint Images 
of the Radio Sky at Twenty-Centimeters (FIRST) radio survey; this is one of the largest BL Lac 
samples yet assembled, and each object emerges with homogeneous data coverage. Each candidate 
is detected in the radio from FIRST and confirmed in SDSS optical spectroscopy to have: (1) no 
emission feature with measured rest equivalent width larger than 5 A; and (2) no measured Ca II H/K 
depression larger than 40%. We subdivide our sample into 426 higher confidence candidates and 75 
lower confidence candidates. We argue that contamination from other classes of objects that formally 
pass our selection criteria is small, and we identify a few very rare radio AGN with unusual spectra 
that are probably related to broad absorption line quasars. About one-fifth of our sample were known 
BL Lacs prior to the SDSS. A preliminary analysis of the sample generally supports the standard 
beaming paradigm. While we recover sizable numbers of low-energy and intermediate-energy cutoff 
BL Lacs (LBLs and IBLs, respectively), there are indications of a potential bias toward recovering 
high-energy cutoff BL Lacs (HBLs) from SDSS spectroscopy. Such a large sample may eventually 
provide new constraints on BL Lac unification models and their potentially peculiar cosmic evolution; 
in particular, our sample contains a significant number of higher redshift objects, a sub-population 
for which the standard paradigm has yet to be rigorously constrained. 

Subject headings: BL Lacertae objects:general — galaxies: active — quasars: general — surveys 



1. INTRODUCTION 

BL Lacs are a rare subclass of Active Galactic Nu- 
clei (AGN) characterized by nearly featureless optical 
spectra, multiwavelength emission, m arked variability 
i and strong polar i zation (e.g., see 'Blandfor d fc ReedllQTSl: 
■ i Kollgaardl [l99l lUrrv fc Padoroni .1995; iPerlman etahl 
\ l2001h . Their spectral energy distributions (SEDs) tend 
to be dominated by synchrotron radiation, with a com- 
ponent due to inverse Compton scattering at higher fre- 
, quencies. The standard paradigm describes BL Lacs in 
the context of AGN unification as FR I radio galaxies 
with their relativistic jets po inted toward the observer 
; (e.g.. iBlandford fc Rees|[l978D . In this scenario the ob- 
served jets are brightened through beaming effects, mak- 
ing BL Lacs powerful probes of AGN jet physics and 
radio-loud AGN in general. Their rarity and lack of 
strong spectral features historically render many tradi- 
tional quasar selection techniques (e.g., UV-excess or 
objective prism surveys) inefhcient. Still, large BL Lac 
samples are eagerly sought but until relatively recently 
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have remained elusive; given their strong multiwave- 
length continuous emission and rarity, broad-band re- 
covery methods over large regions of the sky tend to be 
the most efficient. 

BL Lac studies are sensitive to selection biases and 
historically have been limited by small number statis- 
tics. For example, a now disfavored dichotomy between 
radio selected BL Lacs (RBLs) and X-ray selected BL 
Lacs (XBLs) was initially suggested from a comparison 
between the small and relatively sha. llow but complete 
1 Jy radio sample (jStickel et al.lll99l[) and an X-ray se- 
lected sample from the Einstein Observatory Extended 
Medium-Sensitivity Survey (EMSS: .Stocke et al., 1991) . 
In addition to XBLs appearing to have less rapid vari- 
ability and smaller nuclear to host galaxy flux ratios, 
RBLs and XBLs were thought to populate distinctly dif- 
ferent regio ns of parameter space (e.g., see Figure 5 of 
iPadovani fc G iommi 1995b) . XBLs are also more weakly 
polarized on average, and the majority of XBLs have 
constant polarization positioi i angles (within 20° ) over 
timesc ales of at least 3 year s (iJannuzi et al.|[l99^ . 

But IPadovani fc Giommil ()1995bD argue against two 
distinct BL Lac populations, and they instead propose 
that the observed differences between RBLs and XBLs 
could be explained via different energy cutoffs in the 
beamed synchrotron radiation component of their SEDs. 
They suggest a more physical naming scheme of low- 
energy cutoff BL Lac (LBL) and high-energy cutoff BL 
Lac (HBL) to replace RBL and XBL. LBLs have energy 
cutoffs in the near- infrared/optical, while HBLs have en- 
ergy cutoffs in the UV/X-ray. Most but not all RBLs 
are LBLs in this scenario, and similarly most XBLs are 
HBLs. Deeper surveys have subsequently found many 
objects, aptly named intermediate- energy cutoff BL Lacs 
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(IBL), th at bridge thi s LBL / HBL gap fe.g.. iKock et al. 
1996; Pc rlman et alJ Il996t iLaurent-Muehleisen et aL 
1999i : iColhnge et aLll2005D . thereby demonstrating the 
apparent dichotomy's emergence as a selection effect. 
However, any model advocating a single BL Lac popu- 
lation is still challenged with explaining the phenomeno- 
logical differences between RBLs and XBLs described in 
the previous paragraph. 

Large and deep samples are clearly essential for 
advancing our kno wl edge on BL Lacs (e .g. , see 
Padovani fc Giommil Il995bl: iPerlman et al] 120011 : 
Chen et al and much progress has been made 

over the past 10-15 years: the number of known 
BL Lacs has ballooned from only a couple hundred 
(jPadovani fc Giommil Il995af) to about 1100 cata- 
logued in iVeron-Cettv fc VeronI ([2006) . Numerous 
selection techniques have been used to assemble siz- 
able BL Lac samples. One multifrequency selection 
approach identifies candidates based on flat radio 
spectra, fo llowed by requisit e optica l spectroscop y (e.g., 
Kiihr fc Schmidt 1990; Stickel et a"n [l991: M arch! et all 
19961 ICaccianiga et al.l 2002bD . while another popular 
technique assembles catalogs via the indentification 
of optical counterparts to radio and /or X-ray error 
circles in various s urveys (e.g ., [Stocke et all 119911: 
Wolter et all |1997|; iBadeet ali 119981: iFischer et all 
1998; Laurcnt-Muchlci sen et al] 11998 *: 'Maccacar o et al. | 
1998i: Xaurcnt-Mu ehleisen et aL| 1 19991: iPerlman et al. | 
19981: [Caccianig a et ahP 119991; "iBrinkmann et aLl 
2000'; "Landt ct al.' '2001': 'Caccianiga ct al.' '2002al; 
Beckmann et al. 2003; Padovani ct al. 2007). While 
these combined radio/X-ray selected samples include 
a significant fraction of IBLs in addition to LBLs and 
HBLs, they are of course sensitive to the varying flux 
limits and biases inherent to each multiwavelength 
survey. 

One practical limitation to recovering a large uni- 
form sample from imaging surveys as described above is 
performing the necessary optical spectroscopic follow-up 
for (easily) thousands of candidates, among which only 
< 10^ are typically actual BL Lacs. It has been noted 
that BL Lacs tend to have characteristic X-ray-to-optical 
and optical-to-radio flux ratios differentiating them (at 
least statisti cally) from other radio/X-ray sources (e.g., 
IStocke et al]|1989lll991l: iNass et ali |1996f). This feature 
has been exploited to increase the efficiency of BL Lac 
search algorithms by targeting radio/X-ray sources with 
broad-band flux ratios indicative of known BL Lacs. No- 
tably, the Einstein Slew Survey (IPerlman et al.l Il996f) 
and the Sedentary Multifrequency Sur vey (jGiommi et al.l 
119991 120051 : iPiranomonte et all l2007f l adopted this ap- 
proach to recover large uniform BL Lac samples (^^60 in 
the case of the Einstein survey and '^150 in the final 
Sedentary sample.) While such broad-band color selec- 
tion approaches can recover LBLs, IBLs and HBLs, it 
has been noted that the broad-band SED of BL Lacs can 
vary largely from object to object, and this selection ap- 
proach is prone to missing sor ne BL Lac subpopulations 
potentially in large numbers (lLaurent-Muehleisenlll997l : 
ILaurent-Muehleisen et al.lll999D . It should be noted how- 
ever that the Einstein Slew Survey produced the first 
large BL Lac sample containing a sizable number of 
LBLs, IBLs and HBLs from a single selection technique; 
the Sedentary sample intentionally only attempted to re- 



cover HBLs, and it nicely complements other surveys 
more sensitive to LB Ls (particularly the Deep X-ray 
Radio Blazar Survey: IPerlman et al.l 119981 : iLandt et al.l 
l200lt iPadovani et al.ll200l 

The recent advent of digital large scale spectroscopic 
optical surveys has eased the difficulties of optical follow- 
up incurred by previous BL Lac surveys. A purely opti- 
cally selected sa mple was assembled from the 2df and 6df 
quasa r surveys (iLondish et al.ll2007f l. and lCoUinge et al] 
(|200E |) additionally recovered hundreds of optically spec- 
troscopically selected BL L ac candidates fro m the Sloan 
Digital Sky Survey (SDSS, l^rkeFaDHQOS). These op- 
tically selected samples do not entirely rule out the exis- 
tence of radio-quiet BL Lacs; however, if such a pop- 
ulati on exists the y must be extremely rare (e.g., see 
,Stocke et al.lll99l . 

The ~10^ currently known BL Lacs were derived from 
many different surveys, each with slightly different se- 
lection criteria and biases. In this paper, we present a 
501 object radio sample with homogeneous data coverage 
selected jointly from the SDSS and from the Faint Im- 
ages of the Radio Sky at Twenty-Centimeters (FIRST, 
iBecker et all [T995h survey. This sample is one of the 
largest BL Lac sa mples yet assembled. We also re- 
cently published in I Anderson et al.l (|2007l ) a 266 object 
X-ray sample selected jointly from SDSS spectroscopy 
and from th e Rosat All Sky Survey (RASS. It^es et all 
119991 120001) ■ Our X-ray sample recovered LBLs, IBLs 
and HBLs, and a post-selection correlation with FIRST 
and the NRAO VLA Sky Survey (NVSS) revealed only 7 
objects currently lacking radio detections (possibly, but 
not certainly, because their radio emission lies below the 
FIRST/NVSS catalog flux hmits.) Both our separate 
SDSS/FIRST and SDSS/RASS samples are large with 
similar systematics given their similar recovery methods 
in SDSS. These samples therefore provide an interest- 
ing opportunity for comparing radio and X-ray selected 
BL Lacs in a uniform manner. The sheer size of our 
SDSS/FIRST sample, when supplemented by a rigor- 
ous understanding of selection biases among the different 
wavebands, may eventually permit application to a range 
of BL Lac topics, from AGN unification, to host galaxy 
properties, to BL Lac evolution. 

Our purpose here is to present our radio selected sam- 
ple with some preliminary analysis. We defer detailed 
investigations and comparisons with other samples to a 
later paper. In !j2]we present our selection criteria. Po- 
tential contaminants to our sample are discussed in fj3l 
and our catalog is given in ^l4| Some preliminary results 
and future applications of our sample are discussed in ij5l 
with a summary in Sj6l Throughout we use Hq — 71 km 
s-i Mpc~\ Qm = 0.27, and Aq = 0.73. 

2. SAMPLE SELECTION 

We select BL Lac candidates from the combination of 
the SDSS Data Release 5 (DR5) spectroscopic database 
(|Adelman-McCarthv et a l."2007') in the optica l, and from 
the FIRST survey catalog (|Becker et al.lll995l ) in the ra- 
dio. The SDSS is a multi-institutional effort to image 
10^ deg^ of the north galactic cap i n 5 optical filters cov- 
ering 3800 to 10,000 A (e.g., see iFukugita et all 119961 : 
iGunn et al.l Il998l ) , with follow-up moderate resolution 
spectroscopy (A/AA-1800) of 10^ galaxies lO'^ QSOs, 
and 10^ unusual stars (e.g., see I York et al.ll2000[ ). Data 
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is taken with a special purpose 2.5-meter telescope lo- 
cated at Apache Point Observatory (see iGunn et"al] 
|2006[) , with astrometric accuracy at the ^100 milli-arcsec 
level at the survey limit of r ~ 22 (jPier et all 120031) 
and typical p hotometric prec ision of approximately 0.02- 
0.03 mag Uvezic et al.l 12004 ): 640 simultaneous spec- 
tra are obtained over a 7 deg^ field with a multi- 
fiber optical spectrograph. Fu r ther t echnical details can 
be found in IStoughton et~an (|200l . The SDSS DR5 
spectroscopic database contains spectra for ^10® ob- 
jects over ~5700_deg^, and the DR5 Quasar Catalog 
(jSchneider et al.ll2007l ) contains 77,429 quasars with reli- 
able spectroscopic redshifts. Our selection process iden- 
tifies a very small fraction of objects (< 0.05%) in the 
DR5 Quasar Catalog as BL Lacs (see Q. 

The FIRST radio survey imaged approximately the 
same lO'' deg^ of the North Galactic Cap in the VLA's B 
configuration (5" resolution) at 1.4 GHz. Three- minute 
snapshots were co-added together, an d sources were ex- 
tracted via specialized algorithms fsee lWhite et al.lll997D 
with a catalog flux limit of 1 mJy and positional uncer- 
tainties better than 1". The FIRST survey is by design 
well suited both in sky coverage and depth for match- 
ing optical counterparts in SDSS, and about one-third of 
FIRST sources have detectable o ptical counterparts in 
the SDSS photometric survey (Kimball & Ivczic 2008). 
There are 37,456 objects in the SDSS DR5 spectroscopic 
database that are flagged as FIRST radio sources. 

The target selection algorithms for follow-up SDSS 
spectroscopy are well suit ed fo r reco vering BL Lacs. 
Please refer to [Anderson et al.l (j2003l ) for a thorough 
discussion. Briefly, one goal of the SDSS is to iden- 
tify optical counterparts to ~10'' RASS X-ray sources. 
At least one source with g,r or i < 20.5 within many 
RASS error circles is targeted for spectroscopy. This 
X-ray counterpart algorithm targets many BL Lacs, as 
sources whose SDSS positions additionally match within 
2" to a FIRST source are given highest priority, fol- 
lowed by objects with quasar-like or otherwise strange 
SDSS colors. The SDSS spectroscopic sample is even 
more advantageous for BL Lac recovery in that BL Lac 
spectra are additionally obtained through other quasar, 
serendipity, or gala xy spectroscopi c target selection algo- 
rithms (see Eisens tein et al.l l2001t iRichards et al.l 120021 : 
IStrauss et al.ii2002i) . 

Although one might initially expect that our emphasis 
on radio detection for inclusion would bias our sample 
toward excess LBLs, we do not believe this is an impor- 
tant concern. It is true that by requiring a detection 
in FIRST we are prone to missing some BL Lacs. For 
example, some regions of the sky might not have data 
coverage in FIRST; and, some BL Lacs with SDSS spec- 
troscopy might have been in a low flux state during the 
FIRST observations and would therefore be missed; also, 
we are not highly sensitive to the radio- weak/quiet tail 
of the BL Lac population. While very interesting, radio- 
weak/quict BL Lacs must be extremely rare if th ey exist 
(e.g., see Stocke et al.lll99d iLondish et al.ll2007| ): omis- 
sion of a few would not affect our sample in a statistical 
sense. However, most BL Lacs at the faint flux limit 
for SDSS spectroscopy {g, r, or i < 20.5) will have ra- 
dio fluxes larger than the 1 mJy flux limit of the FIRST 
survey at 1.4 GHz. Figure [1] shows typical SEDs for 
an LBL, an IBL, and an HBL (SEDs were kindly pro- 
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Fig. 1. — Average SEDs for LBLs {solid line), IBLs {dotted 
line), and HBLs {dashed line) in NicDDola et aL (2006), normalized 
to i = 20.5 mag. Approximate flux limits for FIRST in the radio 
(1 mJy at 1.4 GHz), SDSS spectroscopy in the optical (i < 20.5), 
and RASS i n the X-ra y (~2. 5x10-^3 erg s"! cm-^ from 0.1-2.4 
keV, see lAnders on et al. 2Q03) are shown as squares with arrows 
from left to right, respectively. At the optical spectroscopic flux 
limit, we do not expect to miss many BL Lacs by imposing a radio 
flux limit, unless the SDSS object does not overlap with the FIRST 
survey's sky coverage, or if the object was in a low flux state during 
the FIRST observations. Note that applying the RASS flux limit 
could bias our sample toward HBLs. Although we do not formally 
impose an X-ray flux limit to our sample, the RASS flux limit is 
implicit in the SDSS spectroscopic target algorithms. The RASS 
flux limit therefore could affect the makeup of our sample (see JSJ. 
Our conclusion remains unchanged if we choose to normalize each 
SED to g or r = 20.5. 

vided by Elina Nieppola, see iNieppola et al.l l2006f ) as a 
solid line, dotted line, and dashed line respectively; the 
SEDs are normalized to i = 20.5. Overplotted as black 
squares with arrows are typical flux limits for FIRST, 
SDSS spectroscopy, and RASS from left to right, respec- 
tively. So, while it is likely there are some BL Lacs with 
SDSS spectroscopy not detected in FIRST, we do not 
appear to be significantly biasing our sample toward or 
against one particular BL Lac subclass by applying the 
radio flux limit to our selection criteria. Note, how- 
ever, that ap plying the RASS flux limit (taken from 
[Anderson et all l2007l as 2.5x10"^^ erg s~^ cm^^ from 
0.1-2.4 keV) would bias our sample toward HBLs. While 
we do not explicitly apply this flux limit, it does play a 
role in determining which objects are targeted for SDSS 
spectroscopy. Therefore, objects in BL Lac samples de- 
rived from SDSS spect roscopy could be bia sed toward 
HBLs (see [JS] and also IColhnge et al.ll2005D . The rel- 
atively higher RASS flux limit could also explain why 
we recover almost twice a s many BL Lac candidates 
from FIRST/SDSS as An derson et all ([2007.) did from 
RASS/SDSS over the same area of the sky. 

2.1. Selection Criteria 

Following the process detailed below, we impose the 
following criteria for BL Lac selection: 1) the SDSS posi- 
tion matches within 2" to a FIRST radio source; 2) SDSS 
spectroscopy shows no emission line with a rest-frame 
equivalent width (REW) larger than 5 A (see ij2.4p : and 
3) the Ca II H/K break (if present) has a measured de- 
pression C < 0.40. The Ca II H/K depression parameter 
C describes the fractional change in continuum blueward 
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and redward of 4000 A (e.g., see lLandt eralll2002f) : 

iU.r) + [ (Ar) ) ' 

where {fu,b) and {fu,r) refer to the average fluxes per unit 
frequency just blueward (3750-3950 A) and just redward 
(4050-4250 A) of the H/K break respectively. The quan- 
tities (/A,fc) and (/A,r) similarly refer to the average fluxes 
per unit wavelength. These spectral constraints are com- 
monly used in the recent literature, and they are consis- 
ten t with, albeit more restrictive than, those introduced 
bv iMarcha et al.l (jl996). We also take additional pre- 
cautions to reject objects that might formally pass these 
criteria but are not BL Lacs (e.g., stars, post-starburst 
galaxies, obvious broad absorption line quasars, etc.; see 

We note that our adopted selection criteria are em- 
pirically based and thu s somewhat arbitrary (e.g., see 
iRector et al.l Il999l and ICaccianiga et al.l Il999l for rele- 
vant discussions pertaining to selection issues). For ex- 
ample, there is no definitive Ca II H/K break value 
that serves as an obvious transition between beamed 
and unbeamed radio galaxies. And, we could be miss- 
ing some nearby low-luminosit y BL Lacs, as described 
by the Brow ne-Marcha effect ferowne fc Marcii§l ll993t 
iMarcha fcBrowne,1995 ): weakly beamed BL Lacs could 
be so dominated by starlight from their host galaxies that 
they are misidentified as radio galaxies or as clusters of 
galaxies.^ We discuss potential contamination from radio 
galaxies and clusters of galaxies in sections 13.31 and 13.41 
respectively. Also, BL Lacs with significant host galaxy 
contamination can have REW's that are measured to be 
larger than 5 A when referenced to their starlight con- 
taminated continua instead of only their AGN continua 
(jMarcha et al.|[T996[ ). Our selection approach misses such 
objects. Despite these well-studied issues, our adopted 
spectral criteria are also commonly applied in much of 
the recent literature. 

2.2. Initial Query of the SDSS Database 

The SDSS pipeline's measurements of spectral features 
(e.g., redshifts, line equivalent widths, H/K breaks, etc.) 
are of course only useful if the pipeline estimates the 
correct redshift. The pipeline redshifts are quite reliable 
for objects with strong spectral features (e.g., over 98% 
of the pipeline redshifts were accurate for the objects in 
the DR5 Quasar Catalog), but one should not blindly 
trust the redshifts for such weak-featured objects as BL 
Lacs. We therefore placed lenient restrictions on our ini- 
tial database query to ensure high completeness at the 
cost of an initially large false positive rate. We ran the 
following two-part query (followed by manual inspection 
of aU SDSS spectra, see ^TIM : 

1. We return afl SDSS spectra within 2" of a FIRST 
source for which the pipeline has little confidence 
in the measured redshift. This is achieved by se- 
lecting all spectra (that are not blank sky for cal- 
ibration) that the database flags SPECZSTATUS as 
'notjneasured' or 'failed'. This query returned 
589 spectra. 

® Note that this issue is somewhat alleviated by using C = 0.40 
as the H/K break cutoff instead of the more restrictive value of 
0.25 used in some older samples. 



2. For those objects with more confident pipeline red- 
shifts (i.e., SPECZSTATUS ^ 'notjneasured' or 
'failed'), we return all spectra within 2" to a 
FIRST source with the conservative constraints 
that (i) no emission line (of Lya, C IV, C III], 
Mg II, H/3, [O II], or Ha) has a pipeline mea- 
sured REW > 10 A, and (u) the H/K break (if 
present) has a pipeline measured C < 0.41. This 
query returned 1803 spectra. Not every spectrum 
in this set has a reliable pipeline redshift. However, 
as expected, essentially all spectra lacking reliable 
pipeline redshifts also have very weak spectral fea- 
tures, and our query correctly retains these objects. 
Objects with spectral features too strong to be con- 
sidered BL Lacs have correct pipeline redshifts in 
the vast majority of cases and are therefore rejected 
by the query as intended. 

Out of -37,000 total FIRST matches in the DR5 spec- 
troscopic database, our query yielded 2392 spectra as 
potential BL Lac candidates. We do not consider X-ray 
information for selection. However, if an SDSS source 
matches within 1' to a RASS source, we also keep track 
of its X-ray properties listed in the DR5 database. 

2.2.1. Database Query Completeness and Efficiency 

We run two tests to check that our database query 
returns the vast majority of previously-known SDSS 
BL Lacs. F i rst, w e compare our query results to 
ICollinge et all ("2005) 's 183 probable SDSS optically se- 
lected BL Lac candidates that are ad ditionally flagge d as 
FIRST radio sources in their catalog. ICollinge et al.l sim- 
ilarly require C < 0.40 and REW < 5 A, but their 
exact recovery method differs from ours. 172 of their 
183 FIRST/SDSS objects appear in the DR5 database 
as FIRST sources and pass our signal-to-noise cut (see 
§2.3p . Of these 172 sources, our database query fails 
to recover only a single spectrum. We miss this spec- 
trum because the pipeline measured a REW larger than 
10 A for an assumed Mg II feature that is likely just 
noise. 

Second, we perform an identical query for X-ray 
selection by replacing the 2" match to a FIRST 
source with a 1' positional match to a RASS source. 
We then c ompare results from this X-ray query to 
the lAnderson ct al (2007 ) SDSS/RASS X-Ray selected 
AGN sample. lAnderson et al.l queried the DR5 spectro- 
scopic database for SDSS spectra that matched within 
1' to a RASS source, placing no additional constraints 
on the pipeline measured spectral feature strengths or 
redshift flags. Their manual inspection of the resulting 
15,129 spectra yielded 266 X-ray selected BL Lac candi- 
dates (also defined by REW < 5 A and C < 0.40). Our 
revamped X-ray query (i.e., including cuts on equivalent 
width and H/K break) weeds the -10* SDSS/RASS po- 
sitional coincidences in DR5 to —200 objects; we rec ov- 
ered all but 3 spectra included in the lAnderson et ahl BL 
Lac sample. The missed objects include a borderline BL 
Lac candidate with C — 0.4 exactly (as measured by 
lAnderson et aH the SDSS pipeline measured C > 0.41), 
and two examples where the pipeline overestimated a 
single emission line's REW because the line lies in an 
abnormally noisy spectral regime. 

Although our database query is by design inefficient. 
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we conclude that it is excellent in producing a manage- 
able number of candidates while retaining essentially all 
previously-identified BL Lac objects in the DR5 spectro- 
scopic database. 

2.3. SDSS Spectral Signal-to-Noise Cut 

We apply a signal-to-noise cut to the SPSS spectr a 
of potential BL Lacs, similar to iCollinge et al.) (|2005D . 
We calculate the signal to noise over three separate 
500 A wavelength regions centered on 4750, 6250, and 
7750 A (each wavelength region typically contains ap- 
proximately 450, 350 and 300 pixels respectively); these 
spectral regions are typically high-quality and fall within 
the g, r, and i filters respectively. Over each spectral 
region we calculate the signal (S) and the noise (N) such 

that S/N = (E, fx,^/'Jf) / (E. l/'^?)'^', where /a,, and 
(Ti represent the flux density and the estimated uncer- 
tainty in each spectral element i respectively. We retain 
all objects with S/N > 100 in at least one of these spec- 
tral regions. This constraint roughly corresponds to fiber 
magnitudes g < 20.5, r < 20.3, or z < 19.6 for a typical 
spectroscopic plate, and it truncates our list from 2392 
to 1560 objects. 

2.4. Manual Inspection of SDSS Spectra 

We manually inspect the remaining 1560 spectra, in- 
cluding visual verification that our S/N cut did indeed 
only reject objects with poor signal to noise. First, we 
flag each pipeline measured redshift as either reliable, as 
tentative, as a lower limit, or as bad. We define reliable 
pipeline redshifts as those based on at least two spectral 
features at the same pipeline redshift. Tentative pipeline 
redshifts are defined as those showing either only a sin- 
gle spectral feature, typically an emission line assumed 
to be Mg II, or showing multiple weak yet potentially 
real spectral features at the same pipeline redshift. For 
spectra only showing absorption features typically asso- 
ciated with intergalactic medium absorption lines (i.e., 
doublets like Mg II A2796,2803, Fe II A2374,2382, etc.), 
we assign redshifts equal to the lower limit derived from 
the absorption feature and mark those objects' redshifts 
as lower limits. Note that these objects could have red- 
shifts equal to the limit values. Bad pipeline redshifts 
either show no spectral features, or the pipeline simply 
assigned an incorrect redshift (as judged in the visual 
verification process). For objects with tentative and bad 
pipeline redshifts, we attempt to find a redshift that fits 
the spectrum better than the pipeline redshift. If we find 
a better redshift (i.e., >2 spectral features at the same 
redshift,) then we change that spectrum's redshift sta- 
tus to reliable. The small number of spectra that show 
only a single emission line are assigned tentative redshifts 
assuming the lone emission is Mg II. Remaining objects 
with bad pipeline redshifts are then reclassified as having 
unknown redshifts. 

For spectra with reliable and tentative redshifts we 
measure the strongest emission line's REW with IRAF^, 
retaining all objects for which this line's REW < 5 A 
(assuming tentative redshifts are correct.) A handful of 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



objects with redshift lower limits from intervening ab- 
sorption show a single emission feature; we measure this 
emission feature's REW assuming the limit value is the 
correct redshift. We fit any blended lines with multi- 
ple Gaussian components as needed. We do however 
reject candidates with Ha/[N II] complexes that show 
broad Ha, even if our deblended line measurements show 
components with REW < 5 A, under the assumption 
that we underestimate the Ha line flux and overestimate 
the [N II] line flux. Most objects with tentative red- 
shifts have such weak features that they undeniably pass 
the equivalent width cut, especially since they tend to 
have observed-frame equivalent widths < 5 A. Excep- 
tions might be the small number of spectra that show 
only a single emission line. For these objects, we esti- 
mate rest-frame EW by assuming the lone emission is 
Mg II A2800. Objects with adequate S/N but unknown 
redshifts pass our equivalent width cut almost by defini- 
tion. For objects with reliable, tentative, or lower limit 
redshifts, we remeasure the Ca II H/K break parameter 
C (if present) through an automated process, rejecting 
all objects with measured C > 0.40 (assuming tentative 
redshifts and derived lower limits are correct.) We also 
examine the SDSS finding charts of each surviving object 
to ensure none are obviously blended or contaminated by 
a nearby object. 

Finally, we classify each remaining object as either a 
higher confidence BL Lac candidate ('BL') or as a lower 
confidence BL Lac candidate ('BL?'). We require objects 
classified as 'BL' to have slightly higher quality spectra, 
that is S/N > 125 in at least one of the three spectral re- 
gions described in Objects can be classified as 'BL?' 
for any of the following reasons: (1) the largest calculated 
S/N, as described in g231 has 100 < S/N < 125; (2) the 
spectrum passes all criteria, but the object appears close 
enough to another bright object in its SDSS finding chart 
that flux contamination of the fiber spectrum might be 
a concern; or (3) it is unclear where to define the con- 
tinuum near an emission line, and the line's REW is 
measured as larger or smaller than 5 A depending on the 
continuum choice. We classify 426 objects as 'BL' and 75 
as 'BL?'. Some sample spectra for objects classified as 
'BL' and 'BL?' are shown in Figures [2] and [3] respectively. 

2.5. Luminosity and Broad-Band Spectral Index 
Calculations 

We estimate rest-frame monochromatic luminosities 
per unit frequency in the radio, optical and X-ray at 5 
GHz, 5000 A, and 1 keV, respectively. We adopt a radio 
spectral index a^ = —0.27. the average of the 1 Jy radio 
sample (Stick eTet al.lll99lh: and we assum e ao = 1.5 in 
the optical (typical of ICollinge et al.ll2005f s probable BL 
Lac candidates) and ttx = 1.25 in the X-ray (a value 
approximately intermediate between LBLs and HBLs, 
ISambruna et al.lll99i : ISambrunalll997t ). Afl spectral in- 
dices are defined as ~ i/""". We choose the above 
reference frequencies and spectral indices for consistency 
with the bulk of the literature. 

Radio fiux densities at 1.4 GHz are taken from the inte- 
grated FIRST flux densities in the SDSS database. The 
SDSS filter that includes 5000 A in each object's rest 
frame is used to estimate optical fiux densities. We as- 
sume tentative and lower limit redshifts are correct, and 
we use the r filter for spectra with unknown redshifts 
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Fig. 2. — Four sample spectra of higher confidence BL Lac 
candidates, boxcar smoothed with a width of 9. All 4 objects are 
additionally X-ray emitters. Top left: SDSS J082814.20+415351.9, 
an example candidate with a reliable redshift (z = 0.226) and 
prominent host galaxy contamination (Ca II H/K depression 
C = 0.333). Top right: SDSS J111757.24+535554.9, an ex- 
ample of a nearly featureless object with a lower limit red- 
shift (z > 0.720) derived from absorption doublets. Here, 
absorption lines from Fe II A2344, 2374, 2383, 2587, 2600 and 
Mg II A2796, 2804 are all observed at the redshift quoted above. 
Most objects in our sample with redshifts derived from ab- 
sorption doublets typically only show a single absorption fea- 
ture. Bottom left: SDSS J115034.75-I-415440.1, an example 
of a featureless object with unknown redshift. Bottom right: 
SDSS J120412. 11+114555.4, an example of a candidate with a reli- 
able redshift (z = 0.296) and moderate host galaxy contamination 
(C = 0.131). We are careful not to reject objects based on weak 
features near night sky lines such as 5577 A. 
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Fig. 3. — Four sample spectra of lower confidence BL 
Lac candidates, boxcar smoothed with a width of 9. Top left: 
SDSS J082801. 14+231217.6, classified as 'BL?' due to low signal- 
to-noise. This source has an unknown redshift and does not match 
to an X-ray source. Top right: SDSS J095507.88+355100.8, poten- 
tially shows moderate host galaxy contamination (C = 0.231) at a 
tentative redshift z = 0.557. This object is classified as 'BL?' due 
to low signal-to-noise, and it has an X-ray match in RASS. Bottom 
left: SDSS J132558. 94+411500.3, shows moderate host galaxy con- 
tamination (C = 0.223) at a reliable redshift z = 0.310. This source 
is classified as 'BL?' because a broad H« feature can have REW 
larger or smaller than 5 A depending on the continuum level chosen 
near the feature. This source has an X-ray match in RASS. Bottom 
right: SDSS J165329.09+370511.4, classified as 'BL?' due to low 
signal-to- noise. This source has a tentative redshift (z = 1.860) 
and does not have an X-ray match in RASS. 



(as 5000 A would fall in the r filter at the median red- 
shift of o ur sample.) Extinct ion corrected magnitudes 
(using the lSchlegel et al.lll998l dust maps) are converted 
to flux densities using a zero-point of 3631 Jy and effec- 
tive wavelengths 4686, 6165, 7481, and 8931 A for g, r, 
i, and z respectively. We use point spread function (psf) 
magnitudes to reduce potential host galaxy contamina- 
tion in our optical luminosity calculations, but we note 
that we might still be overestimating optical luminosities 
for some nuclei. 

For the X-ray, we use the lStark etaH ([l99l HI maps 
and the colden tool in CIAO (jFruscione et al.l l2006( l 
to estimate the hydrogen column density along each 
sightline. We then use these column densities and the 
Portable, Interactive Multi-Mission Simulator (PIMMS, 
Muka i 1993J to convert ROSAT PSPC X-ray count rates 
to absorbed and unabsorbed broad-band (0.1-2.4 keV) 
fluxes. For each spectrum lacking an X-ray detection we 
estimate its RASS exposure time as the exposure time for 
the X-ray source in the RASS catalog closest to the SDSS 
source's position. Upper limits on RASS count rates are 
then set as 6 counts ( the limit for inclu sion in the RASS 
faint source catalog, iVoges et alj 120001) divided by the 
exposure time, and limits on X-ray fluxes are calculated 
as described above. 

Monochromatic luminosities are then calculated from 
the above radio and optical flux densities and the un- 
absorbed X-ray fluxes, assigning a redshift z=0.297 (the 
median reliable redshift of our 426 higher confidence can- 
didates) to objects with unknown redshifts. We then 
calculate the broad-band spectral indices Ofro, otoxi and 
Oirx-^^ Throughout, we make no attempt to correct lu- 
minosities for the effects of relativistic beaming. 

3. POSSIBLE SOURCES OF CONTAMINANTS 

3.1. Positional Coincidences 

We have very little contamination from random mis- 
matches between SDSS and FIRST source positions. Fig- 
ure [4] (top) shows differences between SDS S and FIRST 
positions for our 501 BL Lac candidates. iBecker et al.l 
(1995) find FIRST positions to be accurate to < 1" at 
the 90% confidence level, and 95% of our SDSS/FIRST 
BL Lac candidates have SDSS source positions within 1" 
of the FIRST radio position. While X-ray emission is not 
a criterion for inclusion in our sample, the bottom panel 
of Figure H pleasantly confirms that we have a low false 
X-ray match rate as well: of our 230 X-ray matches, 89% 
have positional differences less than 25" . This is consis- 
tent with expectations from lVoges et al.l ([1999), who find 
90% of objects in the RASS catalog with stellar optical 
counterparts to have X-ray source positions within 25" 
of the optical counterparts' positions. 

As another test on the level of contamination, we recor- 
relate all 501 BL Lac candidates to the FIRST and RASS 
catalogs after randomly offsetting the SDSS source posi- 
tions by up to ±20" and ±10' respectively (i.e., 10 times 
the corresponding match radii.) We perform the recor- 
relation 10 times, requiring each offset to be larger than 



The broad-band spectral index, for U2 
fined as a^^^a = -"^osiL^^I ^vi)/ 



> 1^1, is de- 
Here, aro = 

log{Lo/ir)/5.08, ciox = —log(Lj:/Lo)/2.60, and Orx = 
— log{Lx/Lr)/7.G8, where Lr, Lo, and Lx are the specific lumi- 
nosities (per unit frequency) at 5 GHz, 5000 A, and 1 keV respec- 
tively. 
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Fig. 4.— Top panel: differences between SDSS and FIRST 
source positions for all 501 BL Lac candidates; 95% of objects 
have SDSS postitions within 1" of each FIRST position. Bottom 
panel: differences between SDSS and RASS source positions for the 
230 candidates with X-ray detections; 89% of objects with X-ray 
detections have SDSS positions within 25" of each R ASS position. 
These percentages are consistent with little contamination from 
positional mismatches between FIRST, SDSS, and RASS. 

each match radius. On average we obtain only 6 ran- 
dom matches to a FIRST or to a RASS source in each 
offsetting test, reaffirming a very low level of potential 
contamination ('^ 1%) due to positional mismatches. 

3.2. Radio Stars 

Since our database queries do not include cuts on ab- 
sorption features other than Ca II H/K break strength, 
we need to remove potential stellar contaminants dur- 
ing manual inspection of the spectra. Manual inspec- 
tion revealed hundreds of radio stars of various spec- 
tral classes identified via common stellar features con- 
sistent with zero redshift. Most of these objects would 
formally fulfill our BL Lac selection criteria, but they are 
rejected since they clearly possess stellar spectra. We de- 
fer substantive discussion of the radio emitting stars for 
a later publication. As an example however, we show 
the spectrum of a previousl y known cataclysm ic variable 
(FIRST J102347.6+003841,lBonditin[lQ0l) in the top 
left panel of Figure [5] that was rejected during manual 
inspection based on emission line strength. 

To test the potential level of stellar contamination in 
our final sample, we examine proper motion informa- 
tion from th e SDSS-I-USNO-B proper motion catalog 
unn et al] [2004) . Proper motions are available for 
~60% of our sample, and we assume objects with proper 



TABLE 1 

Objects with /j > 10.6 milli-arcsec yr^^ and Unknown Redshifts 



SDSS Name 


Class. 




r 


Radio Flux 


(J2000) 


(marcsGc/yr) 


(mag) 


(mJy) 


075547.72+432744.7 


BL? 


11.7 


19.76 


2.63 


094432.33+573536.1 


BL 


52.5 


19.87 


1.86 


094441.47+555752.9 


BL? 


54.5 


19.89 


24.67 


100040.67+531911.7 


BL? 


13.9 


20.36 


2.02 


103845.88+392735.0 


BL 


12.1 


19.23 


7.09 


113115.50+023450.2 


BL 


22.4 


18.53 


13.29 


120355.34+581945.5 


BL 


14.2 


18.57 


41.81 


125032.58+021632.1 


BL 


14.1 


19.21 


468.49 


130146.31+441619.2 


BL 


38.8 


18.64 


53.64 


135158.20+554210.8 


BL 


40.6 


19.22 


83.40 


141208.23+383521.6 


BL? 


13.3 


20.16 


3.35 


164301.06+322104.0 


BL 


25.1 


18.51 


22.21 


165542.81+324420.0 


BL? 


11.2 


18.68 


31.63 



motion information are rep resentative of the rest of our 
sample. iMunn et al.l (|2004f ) find 95% of a sub-sample of 
spectroscopically confirmed quasars with r < 20 in the 
SDSS-I-USNO-B catalog to have proper motions fj, < 10.6 
milli-arcsec year~ ^ . Figure [S] shows a histogram of proper 
motions for 279 BL Lac candidates with proper motion 
detections and r < 20; 91 ± 6% show fi < 10.6 milli- 
arcsec year~^. About half of our BL Lac candidates with 
measured /i > 10.6 milli-arcsec year~^ have reliable or 
tentative redshifts. We therefore do not believe galactic 
sources to be important contaminants to our final sam- 
ple. 

We list 13 BL Lac candidates with large cataloged 
proper motions (/i > 10.6 milli-arcsec yr~^) and un- 
known redshifts in Table [1] We give object names in col- 
umn (1), our classifications as a higher confidence ('BL') 
or lower confidence ('BL?') candidate in column (2), the 
measured proper motions (in milli-arcsec yr~^) from the 
SDSS DR5 database in column (3), psf-magnitudcs in 
the r filter in column (4), and integrated radio fluxes 
from FIRST in mJy in column (5). Only one object 
(SDSS J094432.33+573536.1) has an X-ray detection in 
RASS. These sources could potentially be unusually fea- 
tureless radio stars, and they deserve careful follow-up, 
such as, for example, proper motion measurements with 
very long baseline interferometry. We do not reject these 
sources, however, since we do not have proper motion 
information for every SDSS object and the number of 
objects with large proper motion me asurements is consis- 
tent with random expectations from lMunn et al.) (|2004l ). 

3.3. Radio Galaxies 

Radio emitting galaxies are contaminants returned in 
relatively large numbers by our database query. Most 
contaminating radio galaxies are easily identified through 
emission lines with REW > 5 A and are rejected dur- 
ing manual inspection. Some radio galaxies however can 
be more challenging to recognize. For example, some 
objects show spectra with a weak H/K break due to 
an increased blue continuum not from nuclear activity 
but rather from a recent episode of star formation (e.g., 
E+A ellipticals or other post-starburst galaxies.) Post- 
starburst galaxies with weak emission lines would for- 
mally pass our selection criteria. We therefore reject 
post-starburst galaxies if H(5 shows a REW larger than 
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Fig. 5. — Sample SDSS spectra of objects rejected during 
manual inspection, boxcar smoothed with a width of 9. Top left: 
FIRST J102347.6+003841 (SDSS J102347.68+003841.1), a previ- 
ously known cataclysmic variable (CV, see i|3.2l l. This is the only 
CV recovered by our database query, and it was rejected during 
manual inspection based on emission line strength. Top right : 
SDSS J095842. 64+631845.5 {z = 0.243), an E+A galaxy (see ^T3t . 
This object is rejected as a post star forming galaxy based on H5 
absorption. Bottom left: SDSS J123525. 15+351242.6 (z = 1.870). 
Thi s obj ect is rejected as a broad absorption line quasar (BALQSO, 
see 
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Fig. 6. — Proper motions for 279 BL Lac candidates with r < 20; 
these 279 candidates all have proper motion information available 
in the SDSS database, and 91% show /i < 10.6 milli-arcsec year~^, 
consistent with little stellar contamination in our final sample. 
About half of the objects with fi > 10.6 milli-arcsec year"'^ have 
reliable or tentative redshifts. 

5 A in absorption^ ^, and/or if Ca II H- ^He is significantly 
stronger than Ca II K (as suggested bv lStocke et al.lfl991l 
to reduce galaxy contamination in BL Lac samples). A 
sample spectrum of an E+A galaxy that is rejected dur- 
ing manual inspection based on strong US absorption is 
shown in the top right panel of Figure O 

A small fraction of spectra might show small H/K 
breaks that are not weakened by either nuclear activ- 
ity or star formation, but they could rather be passive 
elliptical galaxies with intrinsically small H/K breaks; 
this does not appe ar to be a strong sou rce of contami- 
nat ion however, as iMarcha et al.l (| 19961 ) point out that 
the lDressler fc ShectmanI ()1987f ) sample of over 700 early 

H(5 absorption is a common proxy for post-star bust galaxies 
(e.g.,[Gotoilal][2003). 



type galaxies has a mean H/K break of C = 0.49 and 
less than 5% show C < 0.40. While requiring the more 
stringent (and historical) criterion of C < 0.25 for inclu- 
sion in BL Lac samples would essentially eliminate any 
concern of contamination from passive elliptical ga lax- 
ies (e.g., see lStocke et al.lll99H ). lMarcha et all (| 19961 ) ar- 
gue that doing so misses some weakly beamed BL Lacs. 
There appears to be a smooth transition between BL Lac 
and radio galaxy spectra, and as noted in SJH the exact 
Ca II H/K break cutoff used to define BL Lacs is thus 
somewhat arbitrary. BL Lac selection criteria are com- 
monly extended to include objects with C < 0.40, as this 
H/K break cutoff would include most beamed AGN, but 
it is not large enough to produce significant contamina- 
tion from elliptical galaxies with intrinsically small H/K 
breaks. We therefore follow the bulk of the recent liter- 
ature and use C < 0.40, acknowledging that this might 
allow a small number of contaminating passive ellipti- 
cal galaxies into our sample. Polarization and variability 
follow-on of objects with 0.25 < C < 0.40 would be help- 
ful for confirming their classification as BL Lacs. 

3.4. Clusters of Galaxies 

See iRector et al.l ()1999() for a discussion on recog- 
nizing BL Lacs from clusters of galaxies. We investi- 
gate in a preliminary fashion potential contamination 
from non-BL Lac objects residing in nearby clusters of 
galaxies by correlati ng our BL Lac candidates to the 
iKoester et al] (I2007D SDSS maxBCG cluster catalog. 
The iKoester et al.l catalog contains 13,823 photometri- 
cally selected clusters with velocity dispersions > 400 km 
from z — 0.1 — 0.3 over 7500 deg^. Photometric 
redshifts are given for each cluster, which are estimated 
to be accurate to Az ~ 0.01 from comparison to avail- 
able spectroscopic redshifts. We perform the matching 
with a fixed projected physical linear aperture of ra- 
dius 0.5 Mpc at the photometric redshift of each cluster. 
Among our BL Lacs, 12 higher confidence candidates and 
1 lower confidence candidate match to a maxBCG clus- 
ter. Of these 13 matches, 4 of the BL Lac candidates 
have no reliable redshift, 5 have reliable redshifts within 
Az ^ 0.01 of the cluster's photometric redshift, and 4 
have reliable, tentative or lower limit redshifts larger than 
the cluster's photometric redshift by at least 0.01. Table 

gives the SDSS name for these 13 BL Lac candidates, 
our classification ('BL' or 'BL?'), the angular separation 
between the BL Lac position and the brightest cluster 
galaxy, each BL Lac candidate's spectroscopic redshift, 
and each cluster's photometric redshift. All objects pass 
our selection criteria, even the more stringent constraints 
aimed toward rejecting post star forming radio galaxies 
(see WS.'S^ : only 3 objects have measured Ca II H/K break 
depressions C > 0.25 (the largest is 0.35). Thus, contam- 
ination from post star forming galaxies or radio galaxies 
with weak H/K breaks residing in clusters is likely small, 
and we retain all 13 objects in our final sample. BL Lacs 
ten d to avoid rich cl usters at low redshifts (z < 0.35, 
see lWurtz et al.llT997[ ). and it is therefore reassuring that 
only a very small number of BL Lac candidates were re- 
turned by our cluster correlation. 

We estimate the expected number of chance superpo- 
sitions by randomly offsetting our BL Lac positions by 
±5 — 50' (the largest angular radius used in our original 
correlation was 4.5') and recorrelating to the cluster cata- 
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TABLE 2 
Matches to maxBCG Clusters 



SDSS Name 


Class. 


Sep. 


BL Lac 


Cluster 


(J2000) 




aresee) 




Zphoto 


082814.20+415351.9 


BL 


27 


0.226 


0.224 


091848.58+021321.8 


BL 


82 




0.278 


092642.18+080300.2 


BL? 


186 




0.116 


101504.13+492600.6 


BL 


126 




0.151 


113118.63+580858.8 


BL 


147 


0.360 


0.167 


122008.29+343121.7 


BL 


64 


0.870"' 


0.246 


124700.72+442318.8 


BL 


63 


1.812 


0.197 


134105.10+395945.4 


BL 


128 


0.172 


0.176 


144248.28+120040.2 


BL 


52 


0.163 


0.159 


161541.21+471111.7 


BL 


68 


0.199 


0.205 


161823.58+363201.7 


BL 


173 


0.734'" 


0.159 


162259.26+440142.9 


BL 


210 




0.122 


163726.66+454749.0 


BL 


73 


0.192 


0.203 



^ rcdshift is a lower limit'" tentative redshift 



log 10 times; we find 5 chance superpositions on average, 
and it is reasonable that the 8 BL Lac candidates in Ta- 
ble [2] with unknown redshifts or redshifts very different 
than each cluster's redshift could be chance superposi- 
tions. Some of the 13 matches (especially those with a 
BL Lac redshift similar to a cluster redshift) might be 
examples of clusters of galaxies hosting a BL Lac. 

3.5. z > 1 BL Lac Candidates and Weak-Lined Quasars 

A follow-on pol arization study of a subset of 42 of 
the ICollinge et ah! (2005) optically-selected BL Lac can- 
didates revealed all nine z > 1 objects in that subset 
lack strong optical polarization (e.g., below 3%), thereby 
raising do ubts as to whether they should be classifed as 
BL Lacs ()Smith et all I20Q7D . We note that these re- 
sults might not be entirely relevant to our sample, as 
only two of their nine z > 1 objects with follow-on 
polarization observations were radio sources. However, 
SDSS has revealed a number of similar objects at very 
high redshift, often termed weak-lined quasars, whose 
optical spectra lack strong emission (like BL Lacs), but 
which also display low polarization and a range of ra- 
dio emissi on extend i ng to quite low values (unlik e BL 
Lacs'): seelFan et al.l (119991 ) . I Anderson et all ()2001[ ). and 
IShemmer et al.l ( 20061 ) for discussions on SDSS weak- 
lined quasars. 

Thirty-one of our higher confidence BL Lac candidates 
and 14 of our lower confidence BL Lac candidates have 
z > 1 (assuming tentative redshifts are correct.) Of our 
31 higher confidence BL Lac candidates with reliable or 
tentative redshifts at z > 1, there are 5 that appear rel- 
atively radio weak (a^o < 0.2), suggesting that some of 
our high-0 candidates are indeed atypical of BL Lacs. 
Given our lack of polarization information we retain all 
z > 1 objects. Polarization, variability, and multiwave- 
length investigations would be enlightening in confirming 
their classification. 

3.6. Broad Absorption Line Quasars 

We also recover several objects that formally fulfill our 
BL Lac criteria but are rejected upon manual inspection 
of their spectra as obvious broad absorption line quasars 
(BALQSOs, see bottom left panel of Figure EJ. 

3.7. An Unusual Class of Radio Sources 



iHall et al] (|2002f ) discovered two objects, 
SDSS J010540.75-003313.9 (z=1.179) and SDSS 
J220445.26+003141.9 (z=1.353), that show a very 
strange continuum drop-off blueward of a Mg II 
absorption doublet; both are FIRST radio objects 
and neither are RASS X-ray sources. Both objects 
also show broad Ha in the near- infrared (Hall et al. 
2008, in preparation). iHall et all were aware of only 
2 other objects previously reported in the literature 
with similar spectra, FBQS 1503-1-2330 (z=0.492) and 
FBQS 1055-^ 3124 (z=0.4(34), b oth of which are also 
radio sources (|White et al.l l200(]f ) . Initially, no entirely 
satisfactory explanatio n was found for these unusual 
objects (see iHall et aLl 's §6.1). H owever, on e promising 
possible explanation suggested by 'H all et al.l is that they 
are extraordinarily unus ual BALQS Os. 

Interestingly, both the lHall et all objects. SDSS J0105- 
0033 and SDSS J2204-h0031, were recovered in our 
database query, and they formally satisfy all our BL Lac 
selection criteria. Our database query also revealed two 
additional radio sources whose SDSS spectra appear very 
similar and may be related: SDSS J130941.36+112540.1 
(z=1.362) and SDSS J145045. 56-1-461504.2 (2=1.877). 
Both are again FIRST sources, and again neither new 
object matches to a RASS source. The SDSS spec- 
tra of all four objects are shown in Figure [71 and Ta- 
ble [3] displays some basic properties. Redshifts for 
SDSS J010 5-0033 and SDSS J2204-I-0031 are taken from 
IHall et all , while the redshifts for SDSS J1309-hll25 and 
SDSS J1450-f 4615 are estimates based on Mg II absorp- 
tion. 

Neither of the two new cases, SDSS J1309-hll25 and 
SDSS J1450-I-4615, ultimately survive as BL Lac can- 
didates during manual inspection: they both appear to 
show broad emission features near Al III/C III. We hence 
exclude the two new cases from our final BL La c catalog, 
although we retain the original two lHall et al.l objects as 
lower confidence BL Lac candidates. Whatever their na- 
ture, this set of radio objects are almost certainly AGN 
and very rare. 

4. FINAL SAMPLE 

Table |4] summarizes the number of objects in our fi- 
nal sample, the number with X-ray detections in RASS, 
and some highlighted redshift information. Thus far, 
from this and previous programs, the SDSS spectroscopic 
survey has already provided several hundred new BL 
Lac candidate identifications; about 60% of our sample 
were prev iously identified as B L Lacs by the SDSS (pr i- 
marily in ICollinge et al.l [20051 or [Anderson et all [20071) . 
and ~20% of our sample were originally discovered by 
non-SDSS observations. Only 33 of our BL Lac can- 
didates appear as qua sars in the DR5 Quasar Catalog 
([Schneider et al.ll2OO70 . Although the DR5 Quasar Cat- 
alog does not include BL Lac objects, some overlap is 
expected between our BL Lac sample and the quasar 
catalog, as our BL Lac definition can include some very 
weak emission line objects. Distributions of observed ra- 
dio flux densities at 1.4 GHz (in mJy), observed optical 
magnitudes (sdss «), and estimated absorbed X-ray fluxes 
from 0.1-2.4 keV (in erg s~^ cm~^, see §2.5p for all 501 
BL Lac candidates are shown in Figure [H 

Data are given in Tables [5| and [51 for our 501 BL Lac 
candidates (426 higher confidence and 75 lower confi- 
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TABLE 3 

Active Galactic Nuclei with Unusual Spectra 



SDSS Name 


rcdshift 


u 


9 


r 


i 


z 


Radio Flux 


(J2000) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mJy) 


010540.75-003313.9" 


1.179 


20.16 


19.04 


17.98 


17.73 


17.49 


4.84 


130941.36+112540.1 


1.362 


21.23 


19.75 


18.57 


17.89 


17.82 


0.83 


145045.56+461504.2 


1.877 


21.00 


20.31 


19.54 


18.80 


18.30 


1.51 


220445.26+003141.9'' 


1.353 


19.78 


18.65 


17.35 


16.76 


16.65 


3.03 



" Discovered bv lHall et"all I I2002D 
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Fig. 7. — Spectra of mysterious radio AGN (poss ibly highly 
unusual BALQSOs; [Hall ct al. 2002) as described in WT\ Spec- 
tra are ordered by increasing redshift and boxcar smoothed with 
a width of 9. None match to a RASS X-ray source. Top left: 
SDSS J010540.75-003313.9 (z = 1.179), published in IHall et all 
This object formally passes our selection criteria, and so it is in- 
cluded as a lower confidence BL Lac candidate, though [Hall et aLl 
have suggested it may be a BALQSO of extreme type. Top right: 
SDSS J220445.26+003141.9 (z = 1.353), published in iHall'et^n 
This object again formally passes our selection criteria, and so 
it is included in our final catalog as a lower confidence BL Lac 
candidate. Bottom left: SDSS J130941. 36+112540.1 (z = 1.362). 
This new object appears similar, but was ultimately rejected 
as a BL Lac based on emission line strength. Bottom right: 
SDSS J145045. 56+461504.2 (z = 1.877). This new object was also 
rejected based on emission line strength. 




-13 -12 -11 
log (erg s"^ cm"^) 

Fig. 8. — Flux distributions for all 501 BL Lac candidates in our 
sample. Top left: Log radio flux densities (in mJy) from FIRST at 
1.4 GHz. Top right: Optical magnitudes (sdss i). Bottom left: Es- 
tima ted l og absorbed X-ray fiuxes from 0.1-2.4 keV in erg s~^ cm~^ 
(see i|2.5p . Open histogram: All 501 candidates. Shaded histogram: 
Limits (denoted with arrows) for the 271 objects lacking RASS de- 
tections. 



dence candidates.) Table [5] presents primarily empirical 
information. Column (1) gives the SDSS source name 
{SDSS JHHMMSS.SS±DDMMSS.S), sorted by right as- 
cension. Column (2) gives our classification as a higher 
confidence BL Lac ('BL') or as a lower confidence BL 
Lac ('BL?') candidate. Column (3) notes whether each 
source is identified as a BL Lac in NED^^. 'S' means 
the source is listed as a BL Lac discovered by SDSS, 
'O' means the source is listed as a BL Lac discovered 
by a survey other than SDSS, and 'N' means it is not 
listed as a BL Lac in NED. Columns (4) and (5) give 
the SDSS coordinates RA and Dec (J2000). Columns 
(6)-(10) present u,g,r,i,z psf magnitudes. Column (11) 
gives the estimated redshift of each object (unknown red- 
shifts are left blank), and column (12) presents our red- 
shift flag: 'R' means we believe the redshift to be reli- 
able, 'T' is a tentative redshift, 'L' means the redshift 
is derived from foreground absorbing gas and should be 
interpreted as a lower limit, and 'U' designates an un- 
known redshift. Column (13) gives each object's mea- 
sured H/K break depression C as derived by us (not the 
SDSS pipeline measure.) Column (14) gives the FIRST 
integrat ed radio flux hsted in the SDSS DR5 database 
in mJy (|Becker et al.lll995HWhite et al.lll997f) . Finally, 
additional comments may appear in column (15). We ex- 
plain in column (15) why each lower confidence candidate 
is classified as 'BL?': 'lowSN' signifies the optical spec- 
trum has measured signal-to-noise 100 < S/N < 125; 
'contam?' means the spectrum might be contaminated 
by a nearby source; and 'ew' marks the objects with 
unclear continua, and emission lines can have measured 
REW larger or smaller than 5 A depending on where one 
chooses to measure the continuum near ea ch line. Not all 
SDSS disco vered BL Lac ca n didat es from lCoUinge et all 
(j2005f ) and lAnderson et all ()2007f) are identified as BL 
Lacs in NED. We thus mark such objects in column 
(3) as not being listed in NED as a BL Lac, but we 
identify all obje cts that appear in lCoUinge et al.l and/or 
lAnderson et al.l as 'C05' and 'A07', respectively, in col- 
umn (15). 

Table [6] emphasizes derived multiwavelength informa- 
tion for each object (see ^2.5p . We repeat SDSS source 
names and redshifts in columns (1) and (2) for easy ref- 
erence to Table El The 1.4 GHz FIRST integrated radio 
flux density (in mJy) is repeated in column (3), with 
the corresponding 5 GHz rest-frame monochromatic lu- 
minosity in column (4). Column (5) gives the extinction 
corrected optical magnitude in the SDSS filter whose 

12 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 
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TABLE 4 
Basic Sample Statistics 





Number of BL 


Number of BL? 


Total 


Number of Objects 


426 


75 


501 


Match to RASS Source 


217 


13 


230 


Reliable Redshifts 


166 


29 


195 


Tentative Redshifts 


58 


17 


75 


Lower Limit Redshifts 


32 


5 


37 


Unknown Redshifts 


170 


24 


194 



center is closest to rest-frame 5000 A, taken from the 
r filter for objects with unknown redshift. Column (6) 
gives the estimated monochromatic optical luminosity at 
rest-frame 5000 A. Column (7) is an X-ray match flag, 
where 'Y' means there is a match to a RASS source, 
and 'N' means there is no match. Column (8) gives 
the X-ray count r ate (in counts s~^) from 0.1-2.4 keV 
(jVoges et al.lll999l) . The estimated monochromatic un- 
absorbed X-ray luminosity at rest-frame 1 keV is given in 
column (9). All monochromatic luminosities have units 
erg Hz~^. The X-ray count rates and luminosi- 
ties presented in columns (8) and (9) are upper limits 
if the X-ray match flag in column (7) is set to N. Fi- 
nally, columns (10), (11), and (12) give the logarithmic 
broad-band spectral indices aro, ctox, and arx respec- 
tively. Lower limits to aox and arx are given if the X-ray 
match flag is set to 'N'. 

As discussed in ^^XTl 171 of the lColHnge et all (|2005D 
optically selected probable BL Lac candidates that 
match to a FIRST radio source are recovered by our DR5 
database query and additionally pass our signal-to-noise 
cut. Two of these objects however are not included in 
our final sample. One object is rejected because it shows 
broad Ha in an Ha/[N II] complex with blended REW 
slightly larger than 5 A. The second object is rejected 
because it shows Si IV+0 IV with REW slightly larger 
than 5 A. 

The [Anderson et all (|2007f) X-ray selected sample con- 
tains 244 objects that match to a FIRST source in the 
DR5 spectroscopic database. ( After a cross-cor relation 
to NVSS, only 7 objects in the An derson et al.l sample 
lack radio detections.) Our database query recovers 241 
of these objects, of which 223 additionally pass our S/N 
cut. Only 2 of these 223 objects do not appear in our 
final radio-selected sample. Both these objects were re- 
jected here based on [O II] emission with REW slightly 
larger than 5 A. Thus, our final sample reassuringly con- 
tains essentially all FIRST sources previously identified 
as BL Lac candidates from SDSS spectroscopy. 

5. DISCUSSION 

The following analysis conservatively focuses on the 
426 higher confidence BL Lac candidates. Results re- 
main virtually unchanged when all 501 candidates are 
included. 

Figure [9] shows a histogram for the 256 higher confi- 
dence BL Lac candidates with either reliable, tentative 
or lower limit redshifts. A sizable fraction (30%) have 
redshifts larger than 0.5 (18% with 0.5 < z < 1.0), 
making this a useful sample for testing BL Lac prop- 
erties at z > 0.5. All objects with z < 1.1 have a 
measured Ca II H/K break, but note that we observe 
a sharp decrease at z 0.5 in the distribution of mea- 
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Fig. 9. — Redshift distribution of the 256 higher confidence candi- 
dates with rehable, tentative or lower hmit redshifts. All redshifts 
with z < l.l are derived from host galaxy features. The relative 
number of objects with measured redshifts decreases drastically for 
redshifts larger than z ~ 0.5. However, a sizable fraction (30%) 
still have measured redshifts larger than 0.5, making this a useful 
sample for constraining higher redshift BL Lac properties. 

sured redshifts. It is possible that some objects with un- 
known redshifts are examples of highly beamed nearby 
BL Lacs; however, given the likelihood of detecting host 
galaxy spectral features from lower redshift objects, it is 
probable that a significant fraction of these objects are 
higher redshift between 0.5 ^ z < 2.2. (Spectra with 
z > 2.2 would start to show Lya forests.) Also, highly 
beamed objects require a low probability geometry and 
may therefore be less likely to be recovered at lower red- 
shifts where a smaller volume is probed. 

A color-color diagram for our 426 higher confidence 
BL Lac candidates is shown in Figure [TOl Objects with 
RASS X-ray detections are marked as blue asterisks, and 
those lacking X-ray detections are shown as black cir- 
cles. Note that X-ray detected and non-detected BL 
Lac candidates cover the same range of optical colors. 
That is, the types of BL Lacs that can show X-ray emis- 
sion are not limited to one specific subclass. The dashed 
lines mark the regio n (g—r < 0.35, r—i < 0.13) where 
iCollinge et al.l (|2005f ) found most stellar contaminants in 
their optically selected BL Lac sample. We find only 11 
objects in this "^ri box" ; none show significant proper 
motion and 8 have tentative or reliable extragalactic red- 
shifts. This affords added confidence that any potential 
stellar contamination to our radio sample is extremely 
small, and it supports the notion that requiring multifre- 
quency information improves BL Lac selection efficiency. 
However, we note that multifrequency approaches are 
not as capable of recovering potentially rare (if extant) 
and extremely interesting BL Lac sub-populations such 
as, for example, radio-quiet BL Lacs. 

Spectra showing larger Ca II H/K depressions tend to 
have redder optical colors on average. The red square in 
Figure fTOl marks the median colors of objects whose spec- 
tra show larger Ca II H/K depressions (0.25 < C < 0.4), 
the yellow triangle marks the median colors of spectra 
showing smaller depressions (0 < C < 0.25), and the 
green star marks the median colors of spectra that do not 
show H/K breaks (either because the break is completely 
washed out by the AGN fiux or because the H/K break is 
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TABLE 5 

Observed Parameters of BL Lac Candidates 



SDSS Name 




in 


RA 


Doc 














Red. 




/l.4 GHz 




(J2000) 


Class. 


NED? 


(J2000) 


(J2000) 


u 


a 


r 


i 




Redshift 


Flag 


c 


(mJy) 


Comm. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


000157.23-103117.3 


BL 


N 


0.48850 


-10.52148 


20.66 


19.74 


18.77 


18.26 


17.86 


0.252 


R 


0.367 


32.97 




000257.17-002447.3 


BL? 


N 


0.73824 


-0.41315 


21.10 


20.63 


19.89 


19.22 


18.85 


0.523 


R 


0.289 


159.13 


lowSN 


002142.25-090044.4 


BL 


N 


5.42608 


-9.01234 


20.15 


19.79 


19.32 


18.82 


18.52 


0.648 


T 


0.217 


41.96 


C05 


002200.95+000657.9 


BL 


O 


5.50396 


0.11610 


20.48 


20.04 


19.28 


18.87 


18.55 


0.306 


R 


0.357 


1.73 


C05,A07 


003808.50+001336.5 


BL 


N 


9.53543 


0.22683 


20.15 


19.70 


19.30 


18.96 


18.63 




U 




89.71 


C05 



Note. — The complete version of this table is in the electronic edition of the journal. The printed version contains only a sample. 



TABLE 6 

Derived Parameters of BL Lac Candidates 



SDSS Name /1.4 ghz X-ray X-ray Counts 

(J2000) Redshift (mJy) log L5 ghz "i^o "^5000 A Match Flag per Sec logLi ^ro otox curx 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


000157.23-103117.3 


0.252 


32.97 


31.82 


18.65 


29.29 


N 


0.019 


25.91 


0.50 


1.30 


0.77 


000257.17-002447.3 


0.523 


159.13 


33.15 


19.15 


29.76 


N 


0.015 


26.63 


0.67 


1.20 


0.85 


002142.25-090044.4 


0.648 


41.96 


32.75 


18.46 


30.16 


N 


0.018 


26.94 


0.51 


1.24 


0.76 


002200.95+000657.9 


0.306 


1.73 


30.71 


19.21 


29.27 


Y 


0.097 


26.85 


0.28 


0.93 


0.50 


003808.50+001336.5 




89.71 


32.40 


19.25 


29.23 


N 


0.019 


26.02 


0.63 


1.23 


0.83 



Note. — The complete version of this table is in the electronic edition of the journal. The printed version contains only a sample. 

In the standard BL Lac beaming paradigm, objects with 
their jets pointed less directly at the observer show less 
non-thermal continua (due to less Doppler boosting), and 
therefore they exhibit more prominent host galaxy fea- 
tmes and redder optical colors. The open histogram at 
the bottom of Figure [TT] shows all 192 higher confidence 
candidates for which we have no measured H/K break 
(either because the H/K break is redshifted out of the 
optical spectrum, or because the object has an unknown 
redshift so it is impossible to locate 4000 A rest-frame.) 
The objects without a measured H/K break tend to be 
bluer on average than those with a larger break (e.g., 
C > 0.2). Also note in Figure [TT] that X-ray detected 
objects are shown as filled circles in the top panel and 
as a filled histogram in the bottom panel; we detect X- 
rays from objects at all ranges of g~r and H/K break 
strength. 

The standard beaming paradigm is further supported 
through an anti-correlation between H/K break strength 
and the log specific luminosity, shown in Figure [T^] for 
the X-ray (top left), optical (top right), and radio (bot- 
tom left). We find linear Pearson correlation coefficients 
r = —0.60 and r = —0.78 for the logarithm of radio and 
optical luminosities respectively; both anti-correlations 
are significant at the >99.999% level. That is, BL Lacs 
appear to be more luminous as we view them more di- 
rectly along their jet axes. Even though there is an ob- 
served trend, the large scatter in Figure [12] would seem 
to prohibit the use of H/K break as a sole proxy for 
predicting accurate luminosities. Also, while our sample 
generally supports the beaming scenario, we note that 
it is still possible that a minority of sources might have 
alternative explanations for exhibiting BL Lac type char- 
acteristics. Detailed investigations of such a possible (al- 
beit likely rare) sub-population are beyond the scope of 
this paper. 

We do not show the luminosity distributions in Fig- 
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Fig. 10. — Color-color diagram of our higher confidence BL 
Lacs. Black circles show higher confidence BL Lac candidates that 
do not match to a RASS X-ray source. Blue asterisks show higher 
confidence BL Lac candidates within 1' of a RASS source. We 
recover X-ray detections at all ranges of optical colors. The green 
star depicts the median colors of objects not showing a Ca II H/K 
break in their spectra. The yellow triangle marks the median col- 
ors of objects with small measured H/K breaks C < 0.25. The red 
square indicates the median colors of objects showing larger H/K 
breaks 0.25 < C < 0.4. As expected, objects with more prominent 
host galaxy features show redder optical co lors on average . The 
dashed lines mark the "gri box" identified bv lCoUinge et aU l|2005l ) 
where most stellar contaminants would be expected. Eight of the 
11 objects recovered in this box have measured redshifts, support- 
ing our conclusion of little stellar contamination to our sample. 
Three objects with g — r > 1.8 and/or r — i > 0.9 due to poor flux 
measurements are omitted from this figure. 



redshifted out of the optical spectrum.) This trend can 
be seen more directly by examining H/K break strength 
vs. color for our higher confidence candidates (Figure[TT|) . 
Our optical spectra are the combination of two distinct 
radiation sources: the host galaxy (typically an elliptical) 
and a typically blue power law continuum from the AGN. 
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Fig. 11. — Ca II H/K depression strength C vs. g — r for our 
higher confidence candidates. Top panel: objects (234) with a 
measured H/K break. The filled circles depict objects that match 
to a RASS X-ray source, while the open circles show those with- 
out an X-ray match. Bottom panel: objects without a measured 
H/K break either because 4000 A rest-fame is redshifted out of 
the SDSS spectra or because the objects have unknown redshifts. 
Open histogram: all 192 higher confidence candidates lacking a 
measured H/K break. The shaded histogram shows the 79 objects 
with an X-ray detection in RASS that also lack H/K break mea- 
surements. Objects with more prominent host galaxy contamina- 
tion (i.e., larger C) tend to be redder, as expected in the standard 
beaming paradigm for the BL Lac phenomenon. Also note that 
objects with X-ray detections are recovered at all ranges of C and 
g — r. Two objects with g — r > 1.8 due to poor flux measurements 
are omitted from this plot. 



Fig. 12. — Multiwavelength luminosities vs. Ca II H/K break 
strength for our higher confldence candidates. Filled circles: ob- 
jects that match to a RASS X-ray source. Open circles: objects 
without an X-ray match (with arrows designating X-ray limits 
when appropriate.) Top left: logarithm of speciflc X-ray luminos- 
ity at 1 keV vs. Ca II H/K break strength C. Top right: logarithm 
of speciflc optical luminosity at 5000 A vs. C. Bottom left: loga- 
rithm of speciflc radio luminosity at 5 GHz vs. C. BL Lacs tend 
to appear more luminous at all observed frequencies when viewed 
more directly along their jet axes (i.e., for smaller values of C). 
Bottom right: logarithm of radio to optical speciflc luminosity ra- 
tio vs. C. The large range of radio to optical luminosity ratios does 
not appear to be an orientation effect and rather probably reflects 
AGN jet physics (e.g., variability and/or differing SEDs among BL 
Lac objects.) Note that we obtain X-ray detections at all values of 
H/K break strength and radio and optical luminosities. 



ure [T2] for objects without a measured H/K break, as 
only 22 such objects have measured redshifts. The aver- 
age X-ray, optical and radio luminosities of 10^* °, lO'^^ '', 
and 10^^'^ (erg s~^ Hz~^) for these 22 objects however 
do lie toward toward the bright end of their respective 
luminosity distributions, as expected. We note though 
that all 22 objects are high redshift (z > 1.222) and 
therefore may be observationally biased toward higher 
luminosities. 

The scatter in optical luminosity vs. H/K break is 
smaller than in the radio or X-ray, as is expected for 
two optical measures. BL Lacs can have diverse SEDs, 
as the peak of the synchrotron component can fall within 
a range spanning close to 10 orders of magnitude in fre- 
quency (e.g., sec Nicppola ct al. 2006). For example, two 
BL Lac jets with identical optical brightnesses can have 
very different radio (or X-ray) luminosities depending on 
the specific shape of each object's SED. This plausibly 
introduces excess scatter into the radio and X-ray panels 
in Figure [T^ especially since we use the same spectral in- 
dices for all objects when estimating luminosities. Also, 
given the non-simultaneity of our multifrequency data, 
variability is likely responsible for some of the observed 
scatter. 

The bottom right panel of Figure [T^] shows the log of 
the ratio of radio to optical specific luminosity vs. mea- 
sured H/K break. Here, we observe no obvious linear 
trend (a 28% two-sided probability of randomly finding 
a stronger correlation), but we do see a relatively large 
scatter. The large scatter likely confirms either variabil- 
ity and/or that BL Lacs do not have a universal SED. 
We omit objects lacking H/K break measurements for 
consistency with the other panels. The lack of an ob- 



vious observed trend with H/K break strength suggests 
that the large observed range of radio to optical luminos- 
ity ratios is not a strong orientation effect and perhaps 
reflects AGN jet physics. 

The bottom right panel of Figure [T^ also illustrates 
one reason that multifrequency surveys should be well 
matched in depth. A radio survey that is shallower than 
an optical survey, for example, will be biased toward re- 
covering objects with larger radio to optical luminosity 
ratios (i.e., extreme LBLs). One needs both a deep radio 
and a deep optical survey to recover objects toward the 
bottom of the bottom right panel of Figure [T^] and to 
obtain a sample with luminosity ratio distributions more 
representative of the underlying population. 

Figures [TUlfT^ demonstrate that BL Lacs display con- 
tinuous distributions in their properties, supporting the 
notion that the observed dichotomy between radio and 
X-ray selected BL Lacs in the 1 Jy and EMSS samples 
was indeed plausibly a selection effect induced by small 
sample size and bright survey flux limits. This is fur- 
ther supported when examining the ~ aox plane in 
Figure [T3] (which shows our 426 higher confidence BL 
Lac candidates.) Blue squares represent objects with X- 
ray detections; arrows designate limits on aox for objects 
lacking a RASS detection. Here, large (small) values 
of aro (ctox) correspond to being radio (X-ray) bright. 
As expected, we observe a continuous transition between 
LBL and HBL objects^'^, with plenty of IBL-type ob- 

HBLs are commonly deflned as curx < 0.75 corresponding to 
Fx/Fr > 1 0~^^ with the X-ra y flux in erg s~^ cm~^ and the radio 
flux in Jy llPadovani fc Giom mi 1995b). The dividing line is also 
sometime s quoted as the similar F^/Fr > 10 ~^'^ where both fluxes 
are in Jy l lWurtzjll99llPerlman et al.|[T99g) . 
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jects bridging the transition. Although this preHminary 
analysis seems to support the notion of a single BL Lac 
population, there are still some differences between his- 
torically radio selected and X-ray selected BL Lacs that 
require further attention. For example, XBLs appear to 
have more starlight contamination from their host galax- 
ies, lower polarization levels, and constant polarization 
position angles (see 21 ■ Also, XBLs and RBLs might 
have very different cosmological evolution. 

We see an obvious and expected bias toward RASS de- 
tected sources preferentially filling the HBL parameter 
space in Figure [131 83% of X-ray emitters appear to be 
HBLs, likely due to the relatively high flux limit of the 
RASS survey. Perhaps counter intuitively however, 66% 
of our entire higher confidence radio sample (and 47% of 
objects lacking RASS detections) still have broad-band 
colors indicative of HBLs. This may be an artifact of the 
SDSS spectroscopic target selection algorithms, which 
are biased toward selecting objects near RASS sources 
for spectroscopy (see 52] and Figure [1]). Keep in mind 
however that the above percentages are tentative for sev- 
eral reasons. About half our data points only have limits 
on their X-ray luminosites; the non-simultaneity of our 
multiwavelength data and our monochromatic luminos- 
ity estimates also introduce excess scatter. 

We also note X-ray detections from some LBLs in 
Figure [131 That is, as is similarly observed in Fig- 
ures [imiT^ X-ray emission is detected from objects at 
nearly all ranges of observed parameters (objects with 
aro ^ 0.2 being the notable exception, see below.) While 
it is certainly plausible that some fraction of radio bright 
objects lacking X-ray detections are real examples of ex- 
treme LBLs, the recovery of X-ray sources at all param- 
eters suggests that we do not have evidence for a signif- 
icant population of X-ray quiet BL Lacs. Deeper X-ray 
observations are necessary for testing this assertion, and 
selection effects and survey biases need to be examined 
thoroughly. 

There is a small population of 7 higher confidence can- 
didates with extremely small a^o < 0.2, of which only a 
single (low-redshift) object shows X-ray emission, imply- 
ing that some BL Lacs might have atypical broad-band 
colors. However, as described in fJSUl five of these objects 
have redshifts z > 1, potentially making their classifica- 
tion as BL Lac candidates suspicious. The other 26 ob- 
jects in our higher confidence BL Lac candidate sample 
with z > 1, however, do not show atypical broad band 
colors. We include all z > 1 objects in our sample with- 
out prejudice, as all formally pass our BL Lac criteria; 
follow-on polarization, variability, and multiwavelength 
studies of these higher redshift BL Lac candidates (or 
weak lined quasars) would be useful. 

Looking forward to future applications, the very large 
size of our homogeneous SDSS/FIRST selected sam- 
ple may (when coupled with a detailed understand- 
ing of selection issues) permit a number of subsequent 
BL Lac studies ranging from the potentially unusual 
evolution of the luminosity function of BL Lacs (e.g. 
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Rector et al.ll20od: iRecto r fc Stock7'2001: Giomm i et alj 
2001t iBeckmann et al.r 2003: Padovani ct al. 2003), to 
tests of the unification paradigm. As just one specific 
future application, consider the standard model that in- 
terprets BL Lacs as FR I radio galaxies with their jets 
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Fig. 13. — Broad-band spectral indices aro vs. aox- Radio 
brightness increases toward the top, and X-ray brightness increases 
toward the left. The blue squares show objects that match to a 
RASS X-ray source while the filled circles depict objects without 
an X-ray match (with arrows designating limits) . We again obtain 
X-ray detections at essentially all ranges of parameters (with the 
radio weak cxro tail being the notable exception). We observe a 
continuous distribution of broad-band spectral indices with plenty 
of IBL and LBL-type objects. However, we do see a potential bias 
toward recovering HBLs (cirx < 0.75) from SDSS spectroscopy, 
a perhaps unexpected result for a radio-selected sample. This is 
maybe an artifact from the SDSS spectroscopic target algorithms. 

pointed toward the observer. Ground and space-based 
observations that resolve BL Lac host galaxies tend to 
support this unification scheme at least for the major- 
ity of objects. An HST survey (primarily with WFPC2 
SNAP observations) found strong direct evidence at 

( Scarpa et all 



lower redshifts supporting this scheme 



2000a^ lTfaV et al.|[2000HFalomo et aL|[2000l:IScan3a et al.1 
I2000b[) . However, it is not yet clear whether it holds 
for every single BL Lac. Suggested modifications range 
from some fraction of BL Lacs being hosted by FR II 
radio galaxies (e.g:..IOwen et al.lll996t lWurtz et al.lll997l : 



IRector fc S tockel l200H iHardcastle et al.l 120031) to some 
fraction Jjeing microlensed fe.g.. 'Ostriker fc V iet"rilll99CII : 

:StockeI£ Rector 1997; Rector & Stockc 2001); any req- 
uisite modifications are likely more significant at higher 
redshifts. Moreover, some caution may be appropriate 
when extrapolating the Urry et al. HST survey results 
to higher redshifts, as only 6 host galaxies with z>0.5 
were actually resolved due to limited WFPC2 SNAP ex- 
posure times. Even if one conservatively disregards our 
75 lower confidence BL Lac candidates and all BL Lac 
candidates with z > 1, our sample still contains a siz- 
able number (46) of higher confidence BL Lac candidates 
with confirmed 0.5 < z < 1 suitable for host galaxy 
investigations. If observed with current and upcoming 
HST instruments, our sample could potentially provide 
an important extension to the 0.5 < z < 1 regime of the 
well-established low redshift unification results. 

6. SUMMARY 

We have assembled a 501 object sample of BL Lac 
candidates selected jointly from SDSS DR5 optical spec- 
troscopy and FIRST radio imaging: there are 426 higher 
confidence candidates and 75 lower confidence candi- 
dates. All 501 objects emerge with homogeneous multi- 
wavelength data coverage from FIRST, SDSS and RASS. 
After imposing a signal to noise cut on the SDSS spec- 
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tra, we require all emission features to have REW < 5 
A and, if present, a Ca II H/K depression C < 0.40. 
Special attention is given to remove radio galaxies, radio 
stars, broad absorption line quasars, and other weak- 
featured objects that might otherwise formally pass our 
selection criteria. Of particular interest, in addition to 
recovering 2 unu sual radio AGN reported previously in 
iHall et aLl (120021 ) . we found 2 additional objects possibly 
of the same class. 

We recover essentially all previous SDSS selected 
BL Lac candi d ates that match to FIRST sources in 
ICollinge et all (|2005f l and lAnderson et~an (|2007f l. in- 
creasing confidence in our selection technique. Of our 501 
BL Lac candidates, approximately 80% are SDSS discov- 
erics (including objects first reported by [A nderso n et al.l 
[2003. Colh ngc et al. 20_0l or I Anderson et a l. 200;^ 

Approximately 60% of our 501 BL Lac candidates have 
either a reliable redshift (derived from at least 2 spectral 
features at the same redshift), a tentative redshift (de- 
rived from a single emission line assumed to be Mg II, 
or from multiple weak spectral features at the same red- 
shift), or a redshift lower limit derived from common 
intervening absorption line doublets. Of our 256 higher 
confidence candidates with spectroscopic redshifts mea- 
sured from SDSS spectroscopy, 30% have z > 0.5, and 
18% have 0.5 < z < 1, making this an excellent sample 
for constraining higher redshift BL Lac properties. 

Almost half of our higher confidence candidates match 
within 1' to a RASS X-ray source. We find X-ray detec- 
tions for BL Lacs at practically all ranges of parameter 
space (i.e., optical colors, multi- wavelength luminosities, 
H/K break strengths, etc.), suggesting that there is not a 
significant population of X-ray quiet BL Lacs. Perhaps 
surprisingly, only about one-third of our higher confi- 
dence radio-selected candidates have broad-band spec- 
tral indices indicative of LBLs. This may be an artifact 
of the SDSS spectroscopy targeting algorithms that could 
bias spectroscopic observations in favor of objects near 
X-ray sources. However, we still recover a large number 
of IBL/LBL-type objects; follow-up of this sample will 
therefore probe a wide range of BL Lac physics. 

Our sample supports the general concurrence that BL 
Lacs display a continuous distribution of properties, but 
we stress that the historical differences between XBLs 
and RBLs still need to be accounted for and require 
further attention. A preliminary analysis seems con- 



sistent with the standard beaming model, in that more 
highly beamed (i.e., more luminous) objects show weaker 
host galaxy features. We also see indirect evidence 
for variability and/or diverse SED shapes, and we note 
that multi-wavelength surveys should be well- matched in 
depth to recover unbiased BL Lac samples. 

Our SDSS/FIRST BL Lac sample constitutes one of 
the largest homogeneous samples to date, with the at- 
tractive attribute that all objects emerge with the same 
high-quality radio, optical, and X-ray data. A thorough 
investigation of the selection effects that shape this sam- 
ple's composition may allow follow-on studies to place 
useful constraints on BL Lac properties ranging from evo- 
lution to unification models, especially extending to the 
high-redshift tail of the BL Lac population. 
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